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ICrutcher. Hakobian fc Troland I (|2009f ) used OH Zeeman observations of four 
nearby molecular dark clouds to show that the ratio of mass to magnetic flux was 
smaller in the ~ 0.1 pc cores than in the ~ 1 pc envelopes, in contradiction to the 
prediction of ambipolar diffusion driven core formation. A crucial assumption was 
that the m agnetic field direction is nearly the same in the envelope and core regions of 
each cloud. IMouschovias fc Tassisl (|2009t ) have argued that the data are not consistent 
with this assumption, and presented a new analysis that changes the conclusions of 
the study. Here we show that the data are in fact consistent with the nearly uniform 
field direction assumption; hence, the original study is internally self-consistent and 
the conclusions are valid under the assumptions that were made. We also show that 
the Mouschovias & Tassis model of magnetic fields in cloud envelopes is inconsistent 
with their own analysis of the data. However, the data do not rule out a more complex 
field configuration that future observations may discern. 
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The relative importance of magnetic fields and turbulence 
in driving the formation of molecular cores is a central ques- 
tion in current star formation theory. Measuring magnetic 
field strengths in order to infer ratios of mass to magnetic 
flux (A//$) has been a focus of observational efforts to an- 
swer this question. However, the only technique for directly 
measuring interstellar magnetic field strengths, the Zeeman 
effect, has only yielded results for the line-of-sight compo- 
nent Blos of the magnetic vector B. A statistical analysis 
carried out with several assumptions has been the standard 
analysis technique for Zeeman results, but this only deter- 
mines the mean or median value of |B| for the observed set 
of clouds, which significantly limits tests of the theory. 

To overcome this limitation, 

ICrutcher. Hakobian fc Troland I (|2009l) (hereinafter CHT) 
carried out OH Zeeman observations toward the envelope 
regions surrounding four molecular dark cloud cores, se- 
lected from a survey of 34 cores for having strong Blos, and 
evaluated the ratio 1Z of M/$ in a cloud core to that in the 
envelope. The goal was to test published strong magnetic 
field models that have uniform initial fields. The value of 
this technique is that published models of core formation 
driven by ambipolar diffusion have strong, regular magnetic 
field morphology such that the unknown angle 6 between 
B and the line of sight is approximately the same in core 



and envelope regions, allowing 9 to be eliminated from the 
expression for 1Z. The CHT analysis depended on the as- 
sumptions that the magnetic field direction was essentially 
uniform in each cloud and that the ratio of strengths of the 
quasi-thermal OH lines between envelope and core gave an 
accurate indicator of the ratio of column densities between 
envelopes and cores. The idealized ambipolar diffusion 
theory of core formation requires 1Z to be approximately 
equal to the inverse of the original subcritical Af/$, or 
TZ > 1. CHT were able therefore to directly test this 
prediction. They found that the probability that all four of 
the clouds have 1Z > 1 is 3 x 10 -7 ; the results are therefore 
significantly in contradiction with the hypothesis that all 
four of these cores were formed by ambipolar diffusion. 

IMouschovias fc Tassisl |2009l) (hereinafter MT) have 
strongly criticized the CHT result, arguing both that the 
CHT analysis is internally inconsistent and that a different 
analysis technique that they apply to the CHT data shows 
that 1Z > 1 is consistent with the data. In this letter we show 
that the CHT analysis is internally consistent, and that the 
MT analysis is itself internally inconsistent. 



2 CONSISTENCY OF THE CHT ANALYSIS 

CHT measured Blos at the four cardinal envelope positions 
(labelled n, s, e, w) surrounding each core, but they inferred 
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Table 1. Observational Results. 



Position Blos (mG) (n — s)/cr (n — e)/a (n — w)/cr (s — e)/cr (s — w)/a (e — iu)/a 



L1448n 
L1448s 
L1448e 
L1448w 


—9 ± 13 
+14 ±8 
-11 ±6 
-7 ±7 


1.5 


0.1 


0.3 


2.4 


0.8 


0.4 


Jz>z ± ( n 

B217s 
B217e 
B217w 


1 1. -4- Q 

+9 ± 13 
+5 ±6 
+6 ±8 


1 A 
1 .1 


1 7 




u. O 




n i 

U. 1 


L1544n 
L1544s 
L1544e 
L1544w 


-3 ±4 
+2 ± 10 

-1 ±4 
+22 ±6 


0.5 


0.4 


3.6 


0.3 


1.7 


3.2 


Bin 
Bis 
Ble 
Blw 


-16 ±6 
-10 ±5 
+0 + 7 
-3 + 6 


0.7 


1.6 


1.4 


1.2 


0.9 


0.4 



the mean Blos in each envelope by fitting simultaneously 
over all the Stokes I and V spectra from the envelope. This 
enabled them to synthesize a toroidal telescope beam that 
sampled the envelope while excluding the core. Blos for 
the core was measured separately by fitting spectra from a 
single beam that covered the core area. They then calculated 
1Z from 

~M Tift ' ' 

+ *lenvelope j ^envelope 

The mass (the OH lines are optically thin) and magnetic 
flux are given by 

MocIAV and $ tx (B L os/cos 6), (2) 

where I is the line intensity, AV the line width, and 9 the 
angle between B and the line of sight. 

With the assumption, based on published strong B 
models (see CHT for references), that 8 cor e ~ enve i ope , the 
cos 6 terms in the numerator and denominator of 1Z cancel, 
and one is left with directly observable quantities only. Thus 
1Z and its uncertainty could be evaluated. CHT evaluated 
the uncertainty with a Monte Carlo analysis that utilized 
the known uncertainties in /, AV, and Blos- This analysis 
hinged on one crucial assumption - that the mean magnetic 
field direction in the envelope immediately surrounding a 
core putatively formed by ambipolar diffusion had the same 
(or nearly the same) direction as in the envelope. MT ar- 
gued that the data are demonstrably inconsistent with this 
assumption and that therefore the CHT analysis is invalid. 

Table 1 and Figure 1 show the relevant data. MT claim 
that the Blos shown in Table 1 show clear variation in 
field strength and perhaps even reversals in field direction 
among the four envelope positions for each cloud. They then 
argue that this variation violates the CHT assumption that 
the field direction is essentially uniform in direction for each 
cloud. 

We also give in Table 1, after the north envelope Blos 
entry for each cloud, each of the six possible differences be- 




sigma sigma 



Figure 1. Left: Histogram plot shows the difference data from 
Table 1, (position 1 - position 2)/cr, in 0.5cr bins; plus signs show 
the normal error function for 24 data points. Right: Same as left, 
except the three differences produced by the L1544w position 
have been removed, and the normal error function is for 21 data 
points. 



tween envelope values of Blos, divided by the la uncer- 
tainty in each difference from the individual measurement 
uncertainties. If MT are correct, these data should show a 
scatter significantly greater than that imposed by the mea- 
surement uncertainties. However, the mean of the 24 differ- 
ences is 1.12 + 0.25, where the uncertainty takes into account 
that there are 16 and not 24 independent measurements. The 
mean for the normal error function is 0.80. Hence, the scat- 
ter in the Blos, including the nominal reversals in direction, 
is only marginally larger than the scatter that would result 
from measurement uncertainty. Most of the larger scatter is 
produced by the single position L1544w; excluding this posi- 
tion, the mean of the remaining 21 differences is 0.87 + 0.18, 
in agreement with measurement error being entirely respon- 
sible for the scatter in the 15 measured Blos excluding 
the L1544w position. The differences are also shown graph- 
ically in Figure 1, along with the normal error function for 
comparison. If the scatter in the Blos were due entirely 
to measurement uncertainty, these two plots would agree. 
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These plots confirm the conclusion drawn from the mean of 
the differences. Even with the L1544w position included, the 
measured differences agree fairly well with the normal error 
function, while if Lf544w is excluded, the agreement is excel- 
lent. Hence, while the Lf 544w position may be anomalous - 
as discussed originally by CHT - there is no statistically sig- 
nificant evidence for the MT claim that the CHT assumption 
of a fairly uniform field direction for each cloud is invalid. 
The scatter in the differences is entirely attributable to the 
measurement uncertainties and not to any intrinsic scatter 
in the Blos- Even if the Lf544 cloud were excluded from 
the CHT analysis, the CHT conclusion that these cores were 
not formed by ambipolar diffusion remains valid. 

MT give an example of possible measurements of f /iG 
and 14 fiG, each with uncertainty 0.1 /xG, and note that the 
mean differs from each value by 2 /xG, not the 0.07 ^G given 
by propagation of errors. However, these 100<r and 140<r ex- 
amples are not germane to the CHT case of roughly 1 — 2a 
measurements. Moreover, CHT did not average the four en- 
velope results for each cloud and obtain the uncertainty by 
error propagation; they synthesized a toroidal beam to sam- 
ple the envelopes and obtained the uncertainties directly 
from the single envelope Blos measurement for each cloud. 



3 CONSISTENCY OF THE MT ANALYSIS 

MT argue that an arbitrarily twisted field morphology (see 
the cartoon shown in MT Figure 1) must be included in the 
analysis for 1Z. Although we have argued above that the data 
do not require that such a morphology is present, let us fol- 
low MT and assume that it is. The MT analysis of the CHT 
data should then be consistent with this proposed model of 
the field morphology - that is, that the angle 9 between the 
core and envelope fields are arbitrarily large. But then the 
MT analysis itself is internally inconsistent with this model. 
CHT defined 1Z in order to eliminate the unknown angle 9 
between the field direction and the line of sight; CHT as- 
sumed that between the core and the envelope of a cloud 
those directions are the same (except for minor differences 
in the 9 that do not significantly affect the analysis, see 
discussion in CHT). That assumption allowed the unknown 
angle 9 between the field and the line of sight, which enters 
as cos 9, to drop out of the ratio 1Z (see Equations 1 and 
2). If the 9s for the four envelope positions and the core of 
each cloud vary greatly, as suggested in MT Figure 1, then 
9s do not drop out of 1Z. For a self-consistent MT analysis, 
each of the five different 9s (core and four envelope) would 
have to be explicitly included in the expression for 1Z. How- 
ever, MT do not do so; such an expression would have the 
five unknown 9s and could not be evaluated. Instead, MT 
use our expression for 1Z with the cos 9s missing. MT stated 
that they were only allowing for different magnitudes in B 
over the four envelope positions, not for different directions 
9. But this assumption is completely inconsistent with the 
astrophysical motivation of strongly twisted field lines (MT 
Figure 1 and discussion) that they give for rejecting the CHT 
analysis and substituting their own. MT offer no astrophys- 
ical explanation for fields at the envelope positions varying 
significantly in strength and perhaps even being antiparal- 
lel while the angle 9 remains invariant. The MT analysis is 



therefore not self-consistent, and cannot be used to analyze 
the CHT data. 

Even within the framework of the MT analysis, it ap- 
pears that the uncertainty in the TZs is overestimated. MT 
considered the variation of the measured Blos as one com- 
ponent of the uncertainty, and then added as a second 
component the measurement uncertainties. However, as we 
showed above, the variation in the measured Blos is consis- 
tent with being due to the measurement uncertainties and 
not to a real variation. MT appear to be doubly counting 
the uncertainties. 



4 CONCLUSION 

CHT concluded that their measurements of the ratios of 
Af/"I> between envelopes and cores did not agree with the 
prediction of the ambipolar diffusion model. Here we have 
shown that the CHT analysis is internally self consistent; 
their conclusions are valid within the framework of the as- 
sumptions they made. The validity of the MT paper rests on 
two pillars: (1) that the CHT data analysis procedure is un- 
ambiguously inconsistent with the data itself, and (2) that 
MT have a superior analysis technique. We have demon- 
strated that neither of these pillars of their paper is correct . 
The conclusions of ICrutcher. Hakobian fc Troland I (|2009l ) 
therefore stand - the observed variations of M/$ from en- 
velope to core are not consistent with the prediction of the 
ambipolar diffusion driven theory of star formation. This 
conclusion does not, of course, rule out the possibility that 
there are structures in magnetic field morphology near dark 
cloud cores; higher resolution and higher sensitivity obser- 
vations would be necessary to investigate this possibility. 

The approach of CHT to test the ambipolar diffusion 
driven model of star formation by measuring the change 
in M/$ between envelope and core is a powerful one that 
should be further exploited, since it reduces uncertainties in 
actual values of magnetic field direction and mass estimates 
by taking ratios. Unfortunately, such experiments will re- 
quire very large amounts of telescope time. However, use of 
the eVLA for OH Zeeman mapping and ALMA for CN Zee- 
man mapping may make it possible to extend this technique 
to smaller scales without requiring such large assignments of 
telescope time. 
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